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1. Introduction 

The nature of the finite temperature phase transition, and the value of the transition temperature 
are basic questions in finite temperature QCD which are still missing final answeres. For QCD with 
2+1 flavors one group ^ finds the deconfining transition, and the chiral transition temperatures are 
separated by 20-30 MeV while another group [Q] claims both temperatures to coincide. Moreover, 
the Brookhaven/Bielefeld collaboration gets for transition temperature = 196(3) MeV, which 
is much higher than the transition temperatures found by the Wuppertal group [||] for the decon- 
fining and chiral transitions - = 170(7) MeV, and = 146(5) MeV, respectively. Both groups 
use rooted staggered fermions, but with different levels of improvement. It has been argued [p 
that the discrepancy is largely due to the rather coarse lattices used by the Brookhaven/Bielefeld 
collaboration. Recently, the Brookhaven/Bielefeld collaboration has extended their calculations to 
lattices of temporal extent A'^r = 8 [@1 and found that with decreasing lattice spacing was shifted 
by 5 — 7 MeV towards smaller values. 

The connection between deconfining and chiral transition has been subject of several phe- 
nomenological considerations. Naively, one would expect the temperature of the deconfinement 
transition to lie below that of the chiral transition, if different at all. This turns out to be the case, 
for example, in the Polyakov-loop extended Nambu-Jona-Lasinio model [||]. More likely is that 
both transitions occur at the same temperature, as Polyakov loop and chiral condensate mix at fi- 
nite dynamical quark masses. The consequence would be a simultaneous enhancement of both the 
chiral and Polyakov-loop susceptibilities along the transition line [|6|, ^ |8[ ^]. 

To clarify the issue, independent investigations of the nature of the finite temperature phase 
transition preferably with different type of the lattice fermionic action are needed (for recent related 
works see [11] (domain wall), [12] (twisted mass) and [ [T3| ] (improved Wilson fermions) ). In this 



work we present results of simulations with Nf = 2 dynamical flavors of nonperturbatively 0{a) 
improved Wilson fermions and plaquette gauge action on lattices of temporal extent Nt = 14, 12, 10 



and 8. Our results were partially reported in [10 |4n. 



2. Definitions and simulation parameters 

The fermionic action for each of the two flavors reads 



-!-£v/(x)[/^(x-a/t) [7Ai + l] V{x-afl) (2.1) 



2a ^ 



za ij^y J 
where P^y is the clover-leaf form of the lattice field strength tensor and 



111 ni^ 

anic = , am = (2.2) 



Kc being the critical value of the hopping parameter. 
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csw 


V = N^Nt 




ro/a 


5.20 


2.0171 


16 8 


0.136050(17) 


5.454(58) 


5.20 


2.0171 


24^10 


0.136050(17) 


5.454(58) 


5.25 


1.9603 


16^8 


0.136273(7) 


5.880(26) 


5.25 


1.9603 


24^8 


0.136273(7) 


5.880(26) 


5.25 


1.9603 


32^2 


0.136273(7) 


5.880(26) 


5.25 


1.9603 


40M4 


0.136273(7) 


5.880(26) 


5.29 


1.9192 


24M2 


0.136440(4) 


6.201(25) 



Table 1: Parameters of the simulation. 



The couplings, lattice volumes and lattice spacings covered by our simulations are listed in 
Table |l|. The scale parameters r^/a have been taken from the zero temperature runs of the QCDSF 
collaboration at the corresponding couplings. They refer to the chiral limit K = Kc. We also list the 
critical hopping parameters jq., which we adopted from QCDSF as well. (For recent relevant work 
see [p3|].) The gauge field configurations were generated on the BlueGene/L at KEK, the RSCC 
cluster at RIKEN, the MVS- 100k at the Joint Computer Center (Moscow), on the SKIF-Chebyshev 
at Moscow State University, as well as on the Altix at HLRN. 

Two-flavor QCD is expected to undergo a second order transition at finite temperature in the 
chiral limit and at very small quark masses. In the chiral limit the order parameter is the chiral 
condensate 

^3 

a = -Y,V{^)v{x). (2.3) 

' X 

For heavy quark masses close to the quenched limit, the theory is known to undergo a first order 
phase transition. In that limit the order parameter is the Polyakov loop 

1 1 ^' 

^ = I ^^^^^^ , L(x) = - Tr n ^4 (x) . (2.4) 

The temperature of the chiral transition is, for general m, identified with the peak position of 
the chiral susceptibility 

X,^{a\ = {a^)-{a)\ (2.5) 
while the peak of the Polyakov-loop susceptibility 

XL^NnL'),,{L'),= {{L')-{Lf) (2.6) 

defines the temperature of the deconfining transition. 

It is expected that the two-flavor theory is in the same universality class as the thi^ee-dimensional 



0(4) Heisenberg model [16], with the external magnetic field and the magnetization being identi- 



fied with the bare quark mass m = am and the chiral condensate d = (a), respectively. The critical 



exponents of this model were found to be [ |17[ ] 



P= 0.537(7), ^=0.206(1). (2.7) 
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Let Tc{m) denote the pseudocritical temperature at finite m, wliicli we define to be the temperature 
corresponding to the position of the peak of the chiral susceptibility 

do 

Xo = ^- (2.8) 
dm 

From the scaling relation connecting the chiral condensate, the dynamical quark mass and the 
temperature in the vicinity of the phase transition we then derive 

1 



Tcirn) — Tc{m = 0) oc . 

Assuming 

oc m, 

we thus expect to find 

Tc{m) - Tc{m = 0) oc m]^^'^^^^ 
for a second order transition at m = 0. A first order transition, on the other hand, would give 

Tc{m) - Tc{m = Q) oc m\. 



(2.9) 
(2.10) 
(2.11) 

(2.12) 



The chiral condensate is related to the average plaquette P by means of a Maxwell relation [p^]. 
The chiral condensate and the plaquette can be found from the partial derivatives of the partition 
function Z: 

1 d 



V dl5 

1 A 

V dm 



InZ 
InZ 



dm,. ^ dcsw ft 
-6P + 2^d-2^^5, 



dl5 



dl5 



2d, 



(2.13) 
(2.14) 



where last term in ( |2.13| ) comes from the clover term and is neglected below because it is suppressed 
by two orders of the lattice spacing with respect to the chiral condensate [|l^, 20, 21]. 

The second derivative d^lnZ/dp dm can be expressed in two different orders, which leads to 
relation: 

dP_ 
dm 



1 djhc 


do 


\ dm 




do 


3 dfi 


dm 


= 3 5^ 


a 


dm 



(2.15) 



called the Maxwell relation. It holds for any lattice size and for all values of j8 and m. Then chiral 
condensate susceptibility can be expressed as 

1 dp 
/I dm ' 

where 



(2.16) 



d rhc d m 



-1 



d^ d^ 
is a finite number. 

We also computed the correlator {Lo)c which can be obtained from the derivative of the aver- 
age Polyakov loop with respect to mass: 

d{L) 



{Lo), 



dm 



(2.17) 
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X 




I ' ' ' ' I 




0.1356 0.1356 0.136 0.1362 0.1364 



0.136 0.1361 0.1362 0.1363 



Figure 1: The chiral susceptibility on the 32^ 12 (left) and 40^ 14 lattice (right) at j3 = 5.25 together with a 
Gaussian fit. 





V 


n)Tc{m) 


XL 


Xa 




5.20 


16^8 


0.682(7) 


2.73(6) 


2.78(6) 


2.81(7) 


5.20 


24^ 10 


0.545(6) 


1.59(8) 


1.59(16) 


1.55(14) 


5.25 


24^8 


0.735(3) 


3.18(4) 


3.17(4) 


3.33(7) 


5.25 


32^12 


0.490(2) 


1.00(11) 


1.05(8) 


1.05(7) 


5.25 


40^ 14 


0.420(2) 




0.59(6) 




5.29 


24^12 


0.517(2) 


1.49(8) 


1.40(9) 


1.3(1) 



Table 2: The pseudocritical temperatures and corresponding pion masses obtained from the peak of Xl, X' 
and {La)c on our various lattices. 



3. Transition temperature 

In Figure [I] we show the chiral susceptibility for lattices with Nt = 12 and 14 corresponding to 
our lowest quark masses. 

In Table ^ we show the pseudocritical temperature Tc{m) and the corresponding pseudocritical 
pion masses, m/, obtained from the peak of the Polyakov-loop susceptibility, the chiral suscepti- 
bility and the correlator (2.17) of L and a, respectively. In Fig. |2|the results together with an ex- 



trapolation to the chiral limit are presented. We find that on all lattices the individual pion masses 
rri^^ coincide with each other within the error bars. 

One can see from the Figure §that the temperature Tc{m) shows an almost Unear behavior in 



the pion mass, in accord with the prediction ( [2. 11[ ) of the 0(4) model. We thus may fit the data by 
the ansatz 

Um)=C + D{rom^)'''\ (3.1) 
The result is shown by the solid curve. Setting the scale by the nucleon mass, the QCDSF collabo- 
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Figure 2: The pseudocritical temperature Tc{m) as a function of pion mass, together with a fit to the power 
m]f''' , according to the three-dimensional 0(4) model. 



ration finds ro = 0.467(15) fm. Using this value, we obtain at the physical pion mass 

roT, = 0.408(5) , T, = 172(3) (6) MeV , (3.2) 
where the first error on Tc is statistical, and the second error reflects the uncertainty in setting the 



scale. This result only slightly differes from our result [14] obtained without 40^ 14 lattice. It is 
in good agreement with the deconfining transition temperature found by the Wuppertal group, but 
lies significantly below the result of the Brookhaven/Bielefeld collaboration. 



4. Conclusions 



We have simulated QCD at finite temperature with two dynamical flavors of nonperturbatively 
improved Wilson fermions on lattices as large as Nt = 14 and lattice spacings as low as 0.075 fm. 
The transition temperature has been computed from the Polyakov-loop susceptibility, the chiral 
susceptibility as well as the correlator of Polyakov loop and chiral condensate. All three tempera- 
tures are found to coincide with each other within the error bars. The critical behavior appears to 
be in accord with the predictions of the 0(4) Heisenberg model, at least as far as the quark mass 
dependence of Tc is concerned. 

Let us note that the Maxwell relation used to compute the chiral condensate susceptibility has 
proven to be a powerful tool in unveiUng the phase structure of clover fermions. 
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